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We present a state-of-the-art evaluation of the polarizability corrections—the inelastic nucleon 
corrections—to the hydrogen ground-state hyperfine splitting using analytic fits to the most recent 
data. We find a value A po i = 1.3±0.3 ppm. This is 1-2 ppm smaller than the value of A po i deduced 
using hyperfine splitting data and elastic nucleon corrections obtained from modern form factor fits. 


Hyperfine splitting in the hydrogen ground state is 
measured to 13 significant figures [l|, 


The proton structure correction inferred from atomic 
hyperfine splitting data by removing other corrections is 


E h{s (e~p) = 1 420.405 751 766 7(9) MHz. (1) 


As = Az + Apoi = -38.58(16) ppm. (4) 


Theoretical understanding of hydrogen hyperfine split¬ 
ting is far less accurate, being at about the part-per- 
million (ppm) level. The main theoretical uncertainty lies 
in proton structure corrections, which are not presently 
calculable from fundamental theory. However, proton 
structure corrections can be calculated as functionals of 
quantities measurable in other experiments, specifically 
as integrals Q over proton form factors measured in elas¬ 
tic electron-proton scattering plus integrals 00, Mi 
over structure functions measured in inelastic polarized 
electron-proton scattering. The quality of the data for 
the latter has improved greatly in recent times, especially 
in the lower momentum transfer region which is impor¬ 
tant for proton hyperfine corrections. In this article, we 
present a state-of-the-art evaluation of the polarizability 
corrections using analytic fits to the most recent data, in 
particular using lower momentum transfer data jg, 0), 10] 
from Jefferson Lab. 

The proton structure corrections, As below, can be 
isolated by taking the experimental values for the hyper¬ 
fine splitting (hfs) and subtracting the other corrections, 
which include QED corrections Aqed, recoil corrections 
A p r , as well as some smaller terms due to hadronic 
vacuum polarization A[[ vp , muonic vacuum polarization 
A(( vp , and weak interactions A(/ eak . One has 0IH 0 

E h{s (e~p) = (l + Aqed + A. P R + A£ vp 

+ A£ vp + A£ eak + As)££, (2) 

where the scale is set by the Fermi energy, given by 


The uncertainty comes mainly from A P R ; we used 5.84 ± 
0.15 ppm. The central value [13] uses the CODATA [5] 
charge radius and the dipole magnetic radius. The un¬ 
certainly encompasses a spread due to other choices of 
these radii. A discussion is in 0 Other quantities were 
taken from 0 Eliminating the large QED corrections 
by using both hydrogen and muonium hfs 0 leads to a 
similar result, with an uncertainty limit 0.18 ppm. 

The structure-dependent correction A5 is convention¬ 
ally split into two terms, Az and A po i. The bulk of the 
first term was calculated by Zemach |2j long ago, and the 
modern expression is 

A z = -2am e r z (l + c>™ d ) , (5) 


where the Zemach radius is 



G e (Q 2 ) 


Gm(Q 2 ) 

1 H“ ftp 



( 6 ) 


Ge and Gm are the electric and magnetic form factors 
of the proton, normalized so that G_e(0) = Gm{ 0)/(1 + 
k p ) = 1, k p = ( g p — 2)/2. The radiative correction 
is estimated in 0 and calculated in 0 for form 
factors that are represented by dipole forms: J^ d = 
(a/37r) [2 ln(A 2 /m 2 ) - 4111/420], For A 2 = 0.71 GeV 2 , 
one finds <51f d = 0.0153; for a A 2 corresponding to the rel¬ 
atively large charge radius in 0 , one would instead get 
0.0150, or a 0.01 ppm change in the hyperfine splitting. 

The second term, A po i, involves contributions where 
the proton is excited j3j„ , 4 la 111 Z| • In the limit where the 
electron mass is neglected except for one overall factor, 


Ep — ol psPp / 
37T (? 




+ m e ) 3 


(3) 


for an electron of mass m e bound to a proton of mass 
trip , magnetic moment p p = (g p /2)(e/2m p ), Lande g- 
factor g p , and fine structure constant a. By convention, 
the exact magnetic moment g p is used for the proton, 
but only the lowest order term, the Bohr magneton pg, 
is inserted for the e~. 


. (y,TTl/e , . * \ 

Apoi =-(Ai + A2), 

7 Tg p m p 

[the prefactor is 0.2265 ppm] with 


Ai 

A2 
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Here F 2 is the Pauli form factor of the proton, 

B 1 = dx(3(r)gi(x,Q 2 ), 

Jo 

P%th 

B 2 = dx0 2 (T)g 2 (x,Q 2 ), ( 9 ) 

Jo 

0{t) = ^ (~3r + 2 t 2 + 2(2 - r) \Jt{t + 1)) and 

/? 2 (t) = 1 + 2r - t(t + 1), (10) 

in which r = e 2 /Q 2 , v is the lab-frame energy trans¬ 
fer, Q 2 is the squared 4-momentum transfer, Xth = 
Q 2 /(2771^771^-1-m^ + Q 2 ), rriTT is the pion mass, and gi and 
g 2 are the spin-dependent structure functions, measured 
in doubly polarized electron proton inelastic scattering. 

Ref. was to our knowledge the first to use g\ }2 data 
to obtain results not consistent with zero for the polariz¬ 
ability corrections. However, Ai and A 2 are sensitive to 
the behavior of the structure functions at low Q 2 . Hence 
with new low Q 2 data available &0 , there is a signifi¬ 
cant possibility that the numerical results for A po i could 
require noticeable revision. 

The integrands of Eq. m converge because there are 
factors that cancel the poles at Q 2 = 0. The function f3 2 
has limiting behavior 


fa (r) 


1/4t = Q 2 /Av 2 
1 


r —> 00 , Q 2 —■> 0 , 
r —> 0, Q 2 —> 00 . 


Given that v is never zero for inelastic processes, even 
for Q 2 —► 0 , and noting that the width of the integration 
region for the Bi is proportional to Q 2 for small Q 2 , one 
sees that the integral for A 2 is well behaved for finite g 2 . 

For the integral Ai to be finite, given 


P(t) 


1 - 5/18r 

0 


OO 
> 0 


( 12 ) 


one needs a cancellation that follows from the Gerasimov- 
Drell-Hearn (GDH) Jig, |2Q] sum rule, 

&1T1 2 

jim 0 QT J dxg ± {x,Q 2 ) = -k£; (13) 

k p = ^( 0 ) is the proton anomalous magnetic moment. 

Further regarding A 2 , there is little g 2 data for the 
proton—there is some from SLAC E155 at higher Q 2 |2lJ 
and there is some Jefferson Lab Hall C RSS data at 
Q 2 = 1.3 GeV 2 under analysis [22]. Hence, the g 2 results 
rely on models. However, the g 2 contributions to the po¬ 
larizability corrections are small because the weighting 
factor 02 in the A 2 integral is generally small within the 
integral. The weighting factor 0 by contrast is on average 
close to 1 . 

Our main results for A po i are detailed in Table [T| We 
evaluated A po i using two different fits to g\ and g 2 and 
several different parameterizations of F 2 (Q 2 ). 


The current best and most up-to-date parameteriza¬ 
tion is the one developed by CL AS EG1 9|. This fit 
begins with the most recent published data | 8 j , which has 
Q 2 > 0.15 GeV 2 , and is based on AO 0] and MAID 24j 
parameterizations of the resonances, the E155 fit 25| 
in the deep-inelastic scattering (DIS) region, and the 
Wandzura-Wilczek 0| form g^ W = —51 + g(y)dy/y 
for g 2 in the DIS region. This fit also gives a good ac¬ 
count of the new data 0 , which has Q 2 down to 0.05 
GeV 2 . JThe other structure function fit is that of Simula 
et al. HI, which is based on data available through the 
year 2001. Our results using this fit are shown in the last 
column of Tabled for one F 2 (Q 2 ). 

For F% , we show results using the parameterizations of 
Kelly |2S| and of Ingo Sick |lSj, and include the dipole 
F 2 (Q 2 ) as a benchmark. Although the dipole fit is a 
traditional standard, it does not fit modern data well 
and results obtained from it are not reliable. The Kelly 
parameterization fits form factor data well overall, and 
for Ge elects to fit the polarization transfer results [§§i]. 
(There has been a discrepancy between the Rosenbluth 
and polarization transfer determinations of the proton 
F 2 . Theoretical analysis 0 is resolving this by suggest¬ 
ing that two-photon corrections to the Rosenbluth de¬ 
termination at non-zero Q 2 will give agreement with the 
polarization method.) The continued fraction parame¬ 
terization of 0 concentrates on the low Q 2 data, and is 
valid for 0 < Q 2 < 0.62 GeV 2 . Beyond this, the F 2 (Q 2 ) 
contributions to the integrals are small and we substi¬ 
tuted the dipole form. This procedure was also used by 
Friar and Sick 0 in their analysis of the Zemach radius. 
Substituting the Kelly parameterization instead made no 
difference on a scale set by our uncertainty limits. 

We show our numerical results in Table |H We have 
split the Q 2 integration into segments [0-0.05] GeV 2 
(where no data exist), [0.05-20] GeV 2 , and > 20 GeV 2 , 
and show contributions from these regions separately. We 
also separate, except for the lowest Q 2 region for Ai, the 
contributions from F 2 , gi, and g 2 . 

Three columns in Table^use the CLAS EG1 model Q 
for gi and g 2 . Errors were assigned to the F 2 contribution 
to Ai using the parameter uncertainties that coherently 
gave the largest error. For Q 2 < 0.05 GeV 2 this error 
was added in quadrature with a 10 % overall systematic 
uncertainty as a guess about the absolute accuracy of the 
data at low Q 2 . The contribution from g 1 was given a 
10% error for the deep-inelastic region ( Q 2 > 20 GeV 2 ), 
and a 50% error (0.05 < Q 2 < 20) dominated by the un¬ 
certainties in the preliminary CLAS data near Q 2 = 0.05 
GeV 2 . The errors on A 2 were taken to be a conservative 
100 %, since there are no significant systematic measure¬ 
ments of g 2 at low Q 2 . To combine errors, we grouped 
the errors due to F 2 , to < 71 , and to g 2 into three sets, 
combined the errors within each set coherently, and then 
combined the errors from the three sets in quadrature. 

When doing the Q 2 < 0.05 GeV 2 part of the Ai in- 
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tegral with the EG1 parameterization, we ensured the 
GDH cancellation by doing a Taylor expansion in Q 2 . In 
terms of the moments 

/*Xth 

r i N 2 } (Q 2 )= x N g h2 (x,Q 2 )dx. (14) 

JO 

(which at low Q 2 satisfy ~ ( Q 2 ) N+1 ) one can ex¬ 
pand Eq. © and obtain B\ = r^ — 10m^rj 2 V(9Q 2 ') + .... 
From the generalized forward spin polarizability |32j, 


7 o (Q 2 ) 


16am 2 f Xth „ 

~Q^~ Jo x 



4 m 2 x 2 \ 

~Q~ ff2 J 


dx , 


(15) 


one obtains —> "/ 0 Q 6 /(16am 2 ) as Q 2 —> 0 , and 

from experiment i33| 70 = [— 1.01 ± 0.08 (stat) ± 0.10 
(syst)]xl0 -4 fm 4 . For the lowest moment, we let Tj 0 " 1 = 
— K 2 Q 2 /(8m 2 ) ± CiQ 4 ± ... at low Q 2 , and integrate the 
low Q 2 part of A- ; keeping terms to 0(Q A ). This gives 


Ai[0,Q?] 



K 2 ± 18 TOpCi 


5 m 2 
4 a 


'7o > <9? , (16) 


where rp is the Pauli radius of the proton. The con¬ 
stant in front of the Q 2 term in T^ reflects the GDH 
sum rule. Calculations in chiral perturbation theory libl 
give Ci = 3.89 GeV -4 , whereas a fit of order Q 6 to 
using the latest CLAS data 0 in the range [0.05-0.30] 
GeV 2 yields ci = 2.95 ± O.ll(stat) GeV -4 . We used the 
experimental value 2.95 GeV -4 for Tabled 

Regarding A 2 , if we assume g 2 = g^ n we find that 
T ( 2 N) = — Nr[ N ' > /(N+1). Therefore, if we naively extrap¬ 
olate this relation to low Q 2 , we get AV0,0.05] = —0.40± 
0.05 using 70 from above. From Ref. (32|, keeping terms 
to 0(Q A ), we find that A 2 [0, Q 2 ] = 3m 2 Q?(yn — Sr,T)/2a. 
Using the MAID 7 r-channel estimate |32j <^lt = 1-35 x 
10 -4 fm -4 , we obtain A 2 [0, 0.05] = —1.4. Both of these 
values would lower the final A po i. We quote the CLAS 
EG1 model (-0.24 ± 0.24) in Table I. 

For the low Q 2 contributions in the Simula et al. j2lj 
fit, we straightforwardly integrated the analytic form. 
The numerical differences from the EG1 result arise be¬ 
cause a Taylor expansion of ( Q 2 ) from Ref. [ 27 ) leads 

to a small curvature parameter ci. 

For all Q 2 , the gi and g 2 fits of have available 
error bands for the gi. To evaluate uncertainties in the 
polarizability corrections due to uncertainties in the gi , 
we recalculated the A, using consistently the largest gi 
and smallest gi. This gives an uncertainty estimate for 
A 2 from g 2 of ±0.37 units, as quoted in Tabled The 
uncertainty estimate for Ai due to g\ and F 2 is ±2.48 
units. 

The uncertainty limits involving g 2 on the proton may 
appear remarkable considering the amount of existing 
data. However, one expects that much of g 2 is due to the 


Wilczek-Wandzura term, which is gotten directly from 
g\. This can be verified from the much larger body of 
data for g 2 on the neutron (using polarized 3 He targets). 

Table GH shows A z and the Zemach radius calculated 
from the form factor sets used above. Table HT1 also gives 
As (data) — Az (calc), i.e.. the value of A po i “desired” for 
consistency between measurement and proton structure 
corrections calculated using a given form factor set. 

In summary, we quote our best value as 

A pol = (1.3 ± 0.3) ppm. (17) 

The earlier value of Faustov and Martynenko [ 3 ] was 
(1.4±0.6) ppm. It is remarkable that this value, based on 
few data agrees with the determination using fits to the 
extensive CLAS data set. We thus corroborate extrac¬ 
tions of hadronic quantities from data as done in GlEl 

We should now focus on modern form factor parame- 
terizations represented by ill or [28|, which fit low Q 2 
data well. Accepting these and Table UTI could lead one to 
desire a larger A po i to reconcile theory and experiment 
for hydrogen hfs. The experimentally-determined A po i 
differs from the value deduced from the measured hyper- 
fine splittings by over 4 ( 6 ) standard deviations of A po i 
for Kelly (Sick). 

We cannot, in our opinion, anticipate that new pro¬ 
ton g 1 or g 2 data will change the evaluation of A po i by 
enough to reconcile the proton structure corrections with 
the measured hydrogen hyperfine splittings. For exam¬ 
ple, to get the requisite A po i by changing the curvature 
parameter we called c\ would require making it about 5 
times larger than the value we used. This is essentially 
unthinkable given data already available. A further look 
at elastic form factor fits that could give a more suitable 
Zemach radius rz could well be warranted. 
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term 

Q 2 (GeV 2 ) 

from 

Kelly’s F 2 

CLAS EG1 

Sick’s F 2 

dipole 

Using Simula et al. gi,g 2 fit 
Kelly’s F 2 

Ai 

[0,0.05] 

F 2 and gi 

0.45 ±0.30 

0.49 ± 0.30 

0.60 ±0.28 

-1.78 ±0.6 


[0.05,20] 

f 2 

7.01 ±0.22 

6.86 ±0.27 

7.12 

7.01 ±0.22 



ffi 

-1.10 ±0.55 

-1.10 ±0.55 

-1.10 ±0.55 

-1.78 ± 1.86 


[20, 00 ] 

f 2 

0.00 

0.00 

0.00 

0.00 



ffi 

0.12 ±0.01 

0.12 ±0.01 

0.12 ±0.01 

0.10 ±0.01 

total Ai 



6.48 ± 0.89 

6.38 ± 0.92 

6.74 ±0.84 

3.55 ± 2.48 

A2 

[0,0.05] 

92 

-0.24 ±0.24 

-0.24 ± 0.24 

-0.24 ± 0.24 

-0.72 ±0.14 


[0.05,20] 

92 

-0.33 ±0.33 

-0.33 ± 0.33 

-0.33 ± 0.33 

-1.14 ±0.23 


[20, 00 ] 

92 

0.00 

0.00 

0.00 

0.00 

total A 2 



-0.57 ±0.57 

-0.57 ±0.57 

-0.57 ±0.57 

-1.86 ±0.37 

Ai ± A 2 



5.91 ± 1.06 

5.81 ± 1.08 

6.18 ± 1.02 

1.69 ±2.51 

Apoi 



1.34 ±0.24 ppm 

1.32 ± 0.24 ppm 

1.40 ± 0.23 ppm 

0.38 ±0.57 ppm 


TABLE I: Contributions to A po i using various models. 


Form factor 

rp (fm) 

r z (fm) Az (ppm) 

residual (ppm) 

Kelly [28] 

0.878(15) 

1.069(13) -41.01(49) 

2.43(52) 

Sick [ 12 ] 

0.871(35) 

1.086(12) -41.67(46) 

3.09(49) 

dipole 

0.851 

1.025 -39.32 

0.74 


TABLE II: Pauli and Zemach radii and Az including 
5 T z for several form factors. The column “residual” is 
As (from hfs data) — Az, i.e., the value of A po i that would 
make theory and data consistent. In this paper, we obtained 
A P oi = (1-3 ± 0.3) ppm using inelastic polarized electron- 
proton scattering data. All uncertainties are calculated from 
the fits except for Sick’s rz and A z, which come from [juii]. 
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